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Abstract: Various moments of the hadronic spectral functions have been employed in 
the determination of the strong coupling from tau decays. In this work we study the 
behaviour of their perturbative series under different assumptions for the large-order be- 
haviour of the Adler function, extending previous work on the tau hadronic width. We 
find that the moments can be divided into a small number of classes, whose characteristics 
depend only on generic features of the moment weight function and Adler function series. 
Some moments that are commonly employed in as analyses from r decays should be avoided 
because of their perturbative instability. This conclusion is corroborated by a simplified as 
extraction from individual moments. Furthermore, under reasonable assumptions for the 
higher-order behaviour of the perturbative series, fixed-order perturbation theory (FOPT) 
provides the preferred framework for the renormalization group improvement of all mo- 
ments that show good perturbative behaviour. Finally, we provide further evidence for the 
plausibility of the description of the Adler function in terms of a small number of leading 
renormalon singularities. 
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1 Introduction 

The precise determination of fundamental parameters of the Standard Model (SM) provides 
one of the most important tests of its internal consistency. In the strong sector, the QCD 
coupling as plays a prominent role and much effort has been devoted to its extraction from 
various observables. The determination from r decays is important, since it provides an 
accurate extraction at low energies, close to the limit of validity of perturbative QCD. 
The general framework for the determination of from the ratio 

= ^ [7 ^ ^-hadrons(7)] ^ ^ 

was developed about 20 years ago [2] . Theoretically, can be expressed as a weighted inte- 
gral of the measured hadronic spectral functions that runs over the hadronic invariant mass 
squared s of the hadronic final state from threshold up to m^. The relevant weight function, 
Wr{x), is obtained from the kinematics of the decay. However, the use of QCD at very low 
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energies is impractical. Therefore, one resorts to a finite-energy sum rule (FESR) where 
the theoretical counter-part of Rr is evaluated as a contour integral in the complex-energy 
plane with \s\ = m^. A particularity of this observable is that non-perturbative effects, 
although small, cannot be neglected. The perturbative QCD result must be supplemented 
with power corrections organised in an operator product expansion (OPE). With the data 
from the ALEPH [3-5] and OPAL [6] collaborations, as well as progress on the theory side, 
the precision on as{mr) is impressive: the advocated uncertainties that two decades ago 
were around 11% [2], are now of the order of 2.5% [5] for the most optimistic analysis. 
However, the results obtained by different groups are sometimes barely compatible, which 
suggests that the details of the different analyses need to be scrutinised. 

In the theoretical description two ingredients are needed. The first of them is the 
perturbative QCD contribution, the second are the non-perturbative effects. At present, 
two theoretical obstacles obstruct progress on the theory side. First, the renormalisation 
group improvement of the perturbative series remains controversial — we return to this 
subject below. Second, the treatment of non-perturbative effects (encoded in the OPE) 
in some of the existent analyses was shown to be inconsistent [7]. A possible solution to 
this problem, proposed in [7-9], is the inclusion of the so-called duality violations (DVs) in 
the analysis framework. These are related to the fact that the OPE fails to describe the 
spectral functions near the Minkowski axis, where resonance effects may become important 
and local quark-hadron duality is violated. In the past, the standard assumption was that 
DVs could be disregarded due to the kinematical suppression of contributions from this 
problematic region. Progress in modelling DVs [10-16] made it possible to include them in 
the as analyses without reliance on external input [8, 9], and to test the above assumption. 

In the present work we focus on the perturbative contribution. One of the main sources 
of uncertainty in the theory of hadronic r decays is the renormalisation group (RG) improve- 
ment of the perturbative series. The most widely employed prescriptions are fixed-order 
perturbation theory (FOPT, see for instance ref. [17, 18]) and contour-improved perturba- 
tion theory (CIPT) [19, 20]. Employing these prescriptions at a finite order in perturbation 
theory leads to differing values for ag. The inclusion of the recently computed correc- 
tion [21] to Rr rendered the discrepancy between FOPT and CIPT even more pronounced. 
Since then, several works have dealt with the RG improvement of the series [18, 22-26]. A 
difficulty common to all these works is that conclusions in favour of FOPT or CIPT (or a 
third prescription) depend on implicit or explicit assumptions on the yet unknown higher 
order coefficients of the Adler function. In particular, the aim of ref. [18] was to construct 
a plausible model for the perturbative series in higher orders incorporating only general 
features of the leading renormalon singularities of the Borel-transformed Adler function. 
This should be sufficient to describe the perturbative coefficients at intermediate and high 
orders, augmented by some polynomial terms to take care of the first few coefficients which 
are not yet dominated by (pre-) asymptotic behaviour. After matching of the model to the 
known coefficients of the Adler function in QCD, the main conclusion of ref. [18] was that 
FOPT is to be preferred over CIPT, since at order and in the region of its smallest 
terms, FOPT provides a closer approach to the resummed series than CIPT. 

How general is this conclusion? A short-coming of ref. [18] and other recent works 
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(with the exception of ref. [24]) is that the analysis was done solely for the kinematical 
weight Wr- This is not entirely satisfactory because the ag determinations employ — and 
often require — several different weight functions, if only to extract Og together with the 
non-perturbative condensates and DV parameters from the same self-consistent analysis. 
In fact, any analytic weight function Wi{x) gives rise to a valid FESR, and different Wi{x) 
emphasise different energy regions of the experimental data and different contributions of 
the theoretical description. In the analyses of ag found in the literature, several different 
moments have been used. The enhancement or the suppression of condensates and DV con- 
tributions have been the guiding principle in choosing these moments. Still, little attention 
has been devoted to the moment dependence of the convergence properties of the pertur- 
bative series and we aim to fill this gap here. In the present paper we therefore pursue the 
FOPT/CIPT comparison using the methods of ref. [18], and ask whether the preference for 
one or the other depends on specific features of the weight function, or the assumptions on 
the Adler function coefficients; whether the kinematic weight is special, or all moments are 
alike. 

The outline is as follows. After setting up the necessary notation, we study how the 
convergence properties of the perturbative expansion depend on the choice of the weight 
function that defines the moment. We try, as much as possible, to remain model independent 
by making use not only of the reference model of ref. [18], but also employing an extreme case 
where CIPT is, by construction, preferred over FOPT for the kinematical weight function. 
We show that with respect to convergence properties and the FOPT/CIPT comparison the 
moments can be divided into several classes, whose global features are manifestations of 
simple properties of the weight functions and assumptions on the Adler function series. We 
conclude that certain weight functions should be more suitable for Og analyses than others. 
We then study the robustness of the model proposed in ref. [18] in the light of the criticism 
presented in ref. [24] and provide further plausibility arguments in favour of the adopted 
procedure. Finally, in the last section we study the consistency between the moments by 
performing simplified Og determinations from single-moment fits. 

2 Theoretical framework 

The total decay rate of the r lepton into hadrons, eq. (1.1), can be separated experimentally 
into three components: the vector, Rr.v, and axial-vector, i?r,A, arising from the decays 
into light quarks through the (■ud)-quark current, and contributions with net strangeness, 
Rt,Si from the (ns)-quark current. Hence 



In determinations of a<j, the focus is on the non-strange contributions, because power cor- 
rections are largest in the strange sector, while they are suppressed by the light u- and 
d-quark masses for Rry ^-nd Rt^a- For this reason, in the following we restrict ourselves to 
the two latter channels. 



Rt = Rt,v + Rt,a + Rt,s ■ 
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The ratios Rry/A can be expressed in terms of integrals of the spectral functions 



Rry/A = 127r5Ew|Kdpy ^ ^1 



mi 



r(0) 



. (2.2) 



In the last equation, Sew is an electroweak correction [27-29], and Vud is the quark-mixing 
matrix element [30] . Theoretically, the relevant two-point functions whose spectral functions 
enter eq. (2.2) are 



(2.3) 



with |r2) being the physical vacuum, and the V and A currents are Jy^^{x) = {u'y^{'y5)d)(x). 
These correlators assume the standard decomposition into transversal and longitudinal 
components which was employed in writing eq. (2.2). 

One then makes use of the fact that the exact correlation functions are analytic in the 
complex s-plane except for a cut along the real axis. This property allows one to write 
eq. (2.2) as a counter-clockwise contour integral along the circle \s\ = sq 



Rv/a(^o) 



67ri5Ew|Kd|^ 



ds 
So 



Wi (s) 



n 



(1+0) 
V/A 



2s ^, 



t(o) 



(2.4) 



|s|=so 



In writing the last equation we have performed two generalisations. First, we are using 
a generalised analytic weight function Wi{s), second, the integral is performed up to an 
arbitrary energy sq < m^. In the notation of eq. (2.4), the particular case of eq. (2.2) 



corresponds to -Ryy^(m^) with sq 



Wt{s) 



and 



1 



mt 



2/ s 
1 + 2—^ 
mi 



(2.5) 



For large enough s, the contributions to n'^"^)(s) can be organised in an operator product 
expansion: a series of local gauge-invariant operators of increasing dimensions times the 
appropriate inverse powers of s. In this framework, the purely perturbative part in the chiral 
limit can be associated with the dimension-zero operator, whereas dimension-2 contributions 
arise from the quark mass corrections.^ The first non-trivial operators arise at dimension 
4, namely, the quark and gluon condensates. The OPE is expected to be well behaved 
along the contour |s| = sq (for sq sufficiently large) except close to the positive real axis. 
Therefore, in the general case, iiy^^(so) obtains a contribution from corrections due to 
the break-down of the OPE close to real s > 0. This term is the aforementioned DV 
contribution. Weight functions Wi{s) that contain one or more zeros at s = sq, such as the 
kinematical w-r in eq. (2.5), tend to suppress the contribution of DVs. 



^In the case at hand, namely u and d quarks only, the dimension-2 corrections are tiny. 
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The different components of -Ry//, can be collected in the following expression 



-Ry/^(-So) = SEwlyudl"^ 



D>2 



(2.6) 



In the last equation, 6^':° and are the perturbative terms, ^ of which 6^} contains the 
Os corrections. Since 6^':^ and 6w] do not depend on the flavour, in the chiral limit they are 
the same for vector and axial-vector correlators, and correspond to the perturbative series 
for the correlator 11^^^^ (s). The contributions from the quark masses, as well as that of 

the operators with D > 2, are encoded in the terms ^l^^y/j^, while the DV contributions 
are represented by S^^y/j^- In this work we are interested in the convergence properties of 

the purely perturbative corrections, and thus our focus is on dw^ . 

The correlator n(^+0) is not RG invariant and contains scale- and scheme-dependent 
contributions. However, the Cauchy integral in eq. (2.4) is insensitive to all s-independent 
terms in the correlators. Without loss of generality, one can work with a renormalisation 
invariant quantity, known as the Adler function, and defined through 

I)(i+o)(,) ^ _^^n(i+o)(^). (2.7) 
as 



Using partial integration, and performing the substitution x = s/sq, one can write 

)(!+( 

r ^ ^v--/ - pert 



= - 27ri / ^ W^{x) I?£+°^(sox) , (2.8) 

|x|=l 



where "pert" denotes the perturbative part of the Adler function in the chiral limit and the 
weight function Wi{x) is obtained from Wi{x) by the integral Wi{x) = 2 dzwi{z). 
In full generality, the perturbative Adler function admits the following expansion: 

,T oo n+1 

4ert^(^) = Y2n^T.<Y.^^-,kL'-\ L^log^, (2.9) 

^ n=0 k=l ^ 

with ttfj, = a{iJ?) = Q!s(//)/7r, fj. is the renormalisation scale, and Nc the number of colours. 
Imposing the RG invariance of the above equation, one may consider as independent only 
the coefficients Cn,i- The other coefficients Cn,k, with k = 2,3,...,n -|- 1, can be obtained 
in terms of the Cn,i and /3-function coefficients.^ (Explicit expressions for some of the 
coefficients Cn^k can be found in eq. (2.11) of ref. [18].) At Nc = Nf = 3 the numerical 
values of the known coefficients c„^i are 

co,i = ci,i = 1, C2,i = 1.640, C3,i = 6.371 [31, 32], C4,i = 49.076 [21]. (2.10) 

Fully analytic results for the coefficients can be found in ref. [21]. Based on a geometrical 
growth of the terms in the perturbative expansion of 6w}, in ref. [18] the estimate 

C5,i ^ 283 (2.11) 



^Henceforth, we omit the so dependence in the terms of the r.h.s of eq. (2.6). 

^We follow the convention of ref. [18], i.e. I3{ay) = ^da^/d^ = X^i,— i /^ftfi/l^^- The first coefficient is 



then /3i = llNc/Q - Nf/3. 
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was put forward for the next term in the series. This estimate was corroborated by the 
model introduced in ref. [18] and we will also employ it in our work. 

Inserting the general expansion of the Adler function, eq. (2.9), into the expression for 

eq. (2.8), yields 

^i?? = f;E^^-^^/^^^(-)l°g'='^(^)<- (2-12) 

71 — 1 k — 1 I I i 

Since the Adler function satisfies a homogeneous RG equation, the above expression for 6w} 
is /x-independent. One has the freedom of setting the scale /i in a convenient way. 

The fixed-order prescription corresponds to fi^ = sq. In this case, the coupling is 
calculated at a fixed scale and can be taken outside the integral. However, the logarithms 
remain to be integrated along the contour. The result can be cast into 

oo n 
n=l k=l 

where the integrals are given by 

^r'"' ^ ^ j'^W.{x)\og^{-x). (2.14) 

\x\=l 

For polynomial moments, these integrals can be performed analytically. Explicit expressions 
for the particular case of the kinematic weight function can be found in ref. [18]. 

In contour-improved perturbation theory [19, 20], the logarithms that remain in the FO 
prescription are summed with the choice /x^ = —sqx before calculating the contour integral. 
This procedure implies that the contour integrals have to be performed over the running 
as in the complex plane 

oo 

^cU = E^".i'^n''""(^o), (2.15) 

n=l 

where the integrals, that can only be computed numerically, are given by 

^n''"''(^o) = ^. (f-W,{x)a^{-sox). (2.16) 
ivrz J x 

\x\=l 

GIFT resums the running of the QGD coupling along the contour of integration. Gonse- 
quently, at each order n, only the coefficient Cn,i enters the expression. 

3 Models for the Adler function 

In order to discuss the behaviour of the perturbative expansion of the spectral moments 
and to compare FO to GI perturbation theory, we need an ansatz for the coefficients c„_i of 
the Adler function, which is the dynamical input common to all moments, beyond n = 4. 
In this section we introduce the models for the series that we use later on. A caveat needs 
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to be spelled out at this point. Going beyond the exactly known coefficients Cn,i requires 
assumptions, usually based on some form of regularity of the series, which might simply be 
wrong. The series might have outliers at some order, and we will never know. The best 
we can do is to state the assumptions clearly, to provide supporting arguments where they 
exist, and to explore the consequences. Two diverging assumptions may be distinguished: 

• The high-order coefficients Cn,i beyond n = 4 are not important, and can be neglected. 
The RG improvement still generates a non-trivial series expansion of the spectral 
moments to all orders through the dependence of Cn^k with A; > 1 on the known c„,^i. 
In a sense this is the assumption underlying CIPT, which assumes that the running 
coupling terms are dominant and therefore should be summed. 

• The high-order coefficients c„^i beyond n = A are essential, since they diverge factori- 
ally for large n, thus overcoming the geometric growth of the running coupling terms. 
Since some knowledge exists on the general structure of this divergence (reviewed in 
ref. [33]), this information can and should be included. 

The two main models that we discuss in this section can be viewed as representatives of these 
assumptions. In addition we also briefly review the result in the large-/3o approximation,^ 
which, since it is based on a well-defined formal limit [Nf — )• — oo) of QCD, provides a 
useful toy model to which we shall return in section 5. In the context of tau decays, this 
toy model has been studied in refs. [34, 35]. 

We start by giving a number of definitions and establishing the notation. We define a 
new function D{s) related to the Adler function by 

— D^^+'\s) = l + D{s) = l + 5]r„a,(V^)"+i. (3.1) 

n=0 

The coefficients c„^i of Dy'^'^^ are related to those of D{s) by c„^i = -K^rn-i- The Borel 
transform of the above series is defined by 

B[D]{t) = (3.2) 

One can then define the Borel integral of the series as (a positive) 

oo 

D{a) = j dte-^/'^ B[D]{t) , (3.3) 



which has the same series expansion in a as D{s) has in as{\/s). The last integral, D{a), if 
it exists, gives the Borel sum of the original divergent series eq. (3.1). An important point 
in the case of the Adler function is that B[D]{t) contains singularities on the positive real 
axis which forces one to adopt a procedure to define the integral D{a). The choice of the 
procedure introduces an ambiguity. We discuss this point in more detail below. 



*For historical reasons, we speak about the "large-/3o" approximation, although in the notation employed 
in this work, the leading coefficient of the /3-function is termed Pi. 
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3.1 Large-/3o model 

In the context of the large-/3o approximation, it has been shown by resumming bubble- 
chain diagrams that the Borel-transformed Adler function has infrared (IR) and ultraviolet 
(UV) renormalon poles at positive and negative integer values of the variable u = f3it/{2Tr), 
respectively [36, 37]. (Except at the value n = 1.) The IR renormalon poles are related to 
the power corrections in the OPE, while the leading UV renormalon dictates the large-order 
behaviour of the series. (For a review see ref. [33]). 

The main result of refs. [36, 37] is that in the large-/3o approximation the Borel trans- 
formed Adler function can be written as [37] 

where the constant C is scheme dependent and cancels the scheme dependence of Og in 
eq. (3.3), such that D{s) is scheme independent. (In the MS-scheme C = —5/3.) In large- 
Pq, all the UV renormalon poles atu=— 1,— 2,..., are double poles, as are all the IR poles 
at u = 3,4, .... The only exception is the IR pole at u = 2, which is simple. This follows 
from the fact that the operator agGG has no anomalous dimension in the large- A'^^ limit. 
The absence of an IR renormalon pole at n = 1 stems from the fact that no dimension-2 
operator contributes to the OPE. The coefficients c„^i in this case can be obtained from 
eq. (3.4) by expanding in u and performing the Borel integral term by term. The first 
12 coefficients can be found in table 1 of ref. [18]. An interesting feature of the large- /3o 
result, in apparent coincidence with the full QCD series (2.10), is that the asymptotically 
dominant sign-alternation from the UV pole at u = — 1 is delayed in the conventionally 
adopted MS-scheme. In intermediate orders the series coefficients are governed by the 
fixed-sign contributions from the u = 2 pole, whose residue is a factor e"^*^ larger. 

3.2 Reference model 

In full QCD we do not have the equivalent of eq. (3.4). On the other hand, the structure 
of the OPE and general RG arguments allow one to determine the position and strength 
of the singularities, which evolve from poles into branch cuts [18, 38, 39], though not their 
residues. In general, the IR and UV singularities are described by the following structures 



B[Df]{u) ^ , l + &i(p-n) + 62(p-n)^ + --- 



1+7 



2 



(j)-u) 

(f^ r - - 1 

i?[DUV](^) ^ P i + b^(^p + u) + b2ip + uf + ... , (3.5) 



(p + u) 



where the constants 7, bi, 7, and 6j of a pole at p depend on anomalous dimensions of 
operators in the OPE as well as /3-function coefficients (the explicit expressions are given 
in section 5 of ref. [18]). When performing the integral (3.3) one needs to circumvent the 
IR singularities along the real axis. A prescription to define the integral is needed which 
introduces an ambiguity in the Borel resummed result. The ambiguity is expected to be 
cancelled by exponentially small terms in ag or, due to the running of the coupling, by 
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power corrections. The treatment of the Borel integral in the presence of these singularities 
is discussed in appendix A of [18]. 

The model for the Adler function constructed in ref. [18] (henceforth called reference 
model, or simply RM) is based on the assumption that it makes sense to merge the exactly 
known low-order behaviour to the leading and sub-leading asymptotics generated by singu- 
larities of the Borel transform. Since the known coefficients do not display the asymptotic 
sign alternating pattern, the leading, first UV singularity should be sufficient. On the other 
hand, if the intermediate orders are governed by fixed-sign behaviour, at least the first two 
IR renormalon singularities should be included in the model. Based on these considerations, 
the ansatz reads 

B[D]{u) = B[DY^]iu) + B[Dl^]{u) + B[Df]iu)+4'^ + d^'^u, (3.6) 

where the renormalon singularities are described by the formulae of eq. (3.5). 

This ansatz is then matched to the known coefficients of the Adler function in QCD as 
follows. First, the known coefficients cs^i, C4,i, and the estimated coefficient cs^i, are used to 
fix the residua of the three renormalon singularities. The polynomial terms are then fixed 
in order to reproduce the lowest order coefficients ci^i and ci^2- The resulting parameters 
are given in eq. (6.2) of [18] and the ffist line of table 3 in section 5, and take "reasonable" 
values. One is then in a position to perform the Borel integration in order to ascribe 
a resummed value to the asymptotic Adler function series. The higher-order coefficients 
Cn,i can be derived and the behaviour of FOPT and CIPT series can be compared to the 
resummed one. 

The main conclusion of ref. [18] is that under the above assumptions FOPT is clearly 
preferred over CIPT for Wr- The CIPT series displays a faster convergence but fails to 
give a good approximation to the Borel resummed result in the sense of an asymptotic 
series. From FOPT, on the other hand, it is possible to extract a good approximation to 
the Borel resummed value in spite of the slower convergence of the series. The reason for 
this observation can be traced back to cancellations that are missed by the CIPT series. 
To understand this we rewrite ^ of eq. (2.13) as 

oo 

^fL. = E h-i'^r + at^] «(^o)" , (3.7) 



with 



n=l 



=5:A:c„,,J™r. (3.8) 

k=2 

j;tree 



The Cn,i series is simply the Adler function series multiplied by 6Jj'r^, while the contour 
integration of the running coupling effects is fully contained in the gjj"*^ series. In eq. (3.7) 
the tree-level contribution arises because J^'""' = Wi{0) = FOPT treats the c„,i and 
the g^''^ series on an equal footing. Comparing this decomposition with the CIPT result, 
eq. (2.15), one observes that in CIPT the gli ^ series is resummed to all orders while the Cn i 
series is used only up to a finite order n. An important mo del- independent feature of the 
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Reference model 



Alternative model 



n 






(c„ + .glr^' ) /c„ 






(c„ + .gn""' ) /C„ 


4 


49.1 


78.0 


2.59 


49.1 


78.0 


2.59 


5 


283 


307.8 


2.09 


283 


307.8 


2.09 


6 


3275.4 


-807.3 


0.75 


2148.3 


-807.3 


0.62 


7 


18,758 


-10,398 


0.45 


11,801 


-34, 489 


-1.92 


8 


388, 442 


-329,054 


0.15 


150,508 


-592,196 


-2.93 


9 


919,121 


-232,718 


0.75 


215,264 


-5.1 X 10^ 


-22.8 


10 


8.4 X 10^ 


-7.3 X 10^ 


0.12 


2.4 X 10^ 


-6.4 X 10'^ 


-1.69 



Table 1. Cancellations between the c„ and the gn ^ series for orders 4 < n < 10 for the kinematical 
moment, w-r, in the reference model and in the alternative model, eqs. (3.6) and (3.9). 



QCD series, which follows from the OPE and the form of the moment weight function, is 
that there are large cancellations of n! divergences between the c„^i and g^''^ series. These 
cancellations are particularly strong when the series is dominated by the u = 2 singularity, 
and for the kinematic weight.^ In such a scenario, it is mandatory to combine Cn,i and g^^^ 
order by order in n, lest the cancellations do not take place. Since CIPT treats the orders 
incoherently, it misses the cancellations and runs into the sign alternating asymptotic regime 
earlier. In FOPT, on the other hand, the cancellations suppress the divergence and allow 
FOPT to approach the Borel result. In table 1, we show as an example the cancellations 
in the case of the kinematical moment Wr- However, note that they are not imposed in the 
RM. Rather, the matching procedure to the QCD series gives the expected weight to the 
leading IR pole. If the residue of the IR pole u = 2 had turned out to be tiny, this would 
have made the cancellations almost non-existent. 

3.3 Alternative model 

To make this feature clearer, we can artificially suppress the leading IR pole. Let us 
consider a model for the Borel transformed Adler function where the IR singularity at 
u = 2 is removed and another at u = 4 is added: 

B[D]{u) = B[DY^]iu) + B[Df]{u) + B[Df']iu)+4'^ + d^'^u. (3.9) 

In this model, the aforementioned cancellations do not take place by construction. We refer 
to this model as the alternative model (AM). An analogous matching procedure can be 
carried out yielding the following values for the parameters: 

4^ = 66.18, 4^ = -289.71, = -5.21 x 10-^ 

40 = 2.15, 40^401x10-1. ^"^-^"^ 

Here, as we show in the sequel, CIPT is able to approach the Borel resummed result while 
FOPT exhibits oscillations around this value. Table 1 shows that the cancellations between 

^In large-/3o this is shown analytically for Wr in ref. [18]. An important point discussed in the next 
section is the moment dependence of the cancellations. 
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the Cn and gh, series no longer take place in this model. The table also shows a slower 
growth of the Cn,i in this model, and a dominance of the gn terms up to n = 10, which 
therefore realises a situation where running coupling effects are dominant. We use this 
model as an example where CIPT is, by construction, superior to FOPT, at least for the 
kinematical weight. This provides a way to assess a possible model dependence in our 
conclusions. Nevertheless, we emphasise that we find it unlikely that the Adler function in 
QCD behaves as the AM, since there is no known mechanism that would naturally suppress 
the u = 2 singularity. 

4 Moment analysis 

The determination of as and condensates from the analysis of r hadronic spectral functions 
is based on sum rules obtained by equating eqs. (2.2) and (2.4). In the former, to perform 
the integral along the real axis, the experimental spectral functions are used. In eq. (2.4), 
the theoretical description of the correlators is employed in the contour integration. An 
important aspect of these sum rules is that one still has the freedom of choosing any 
analytic weight function Wi{x), as well as any point sq < (as long as sq is large enough 
for the OPE and the perturbative expansion to make sense). On the experimental side, it 
is obvious that a given weight function enhances the regions of the spectrum where it has 
peaks. On the theory side, the relative contributions of the different 5's in eq. (2.6) are 
strongly dependent on the choice of Wi{x). For example, as already mentioned, moments 
of functions Wi with zeros at s = sq suppress ■ These are known in the literature under 
the name pinched moments. A monomial term of the type in Wi, on the other hand, 
implies the monomial x^^^ in Wi defined after eq. (2.8), which enhances (or, rather, does 
not suppress) the contribution of the condensate of dimension D = 2{k + 1) as well as the 
factorial divergence from the IR renormalon singularity at u = k+1. Let us analyse, as an 
illustrative example, the kinematical moment 

Wr = (1- xf{l + 2x) = l-3x'^ + 2x^. (4.1) 

It receives its larger contributions from 6^'^'^, 6w}, 6w}, and 6w}. The first two arise mainly 
from the 1 in Wr, whereas 6w} and 6w} arise from the terms — and 2x^, respectively. 
The double zero at j; = 1 suppresses , while the absence of other monomial terms 
suppresses the condensates with D = 4 as well as with D > 10. (The mass corrections, 
6w^, are negligible due to the smallness of the quark masses.) 

In order to extract as, a number of condensates, and the DV parameters from the 
data sets one needs more than one observable. It has become standard to use a set of 
several weight functions Wi in order to perform a combined fit to their — not statistically 
independent — moments. In table 2 we collect the weight functions investigated in this 
work. Most of them have been employed in at least one of the recent analyses of hadronic 
r spectral functions. In this table, the first five rows are the building blocks for the other 
polynomial weight functions. The second set are pinched weight-functions that contain a 1 
followed by powers of x. The third block contains pinched weight-functions that do not have 
the 1 and start directly with some power of x. The idea behind the moments that were used 
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Table 2. Weight functions investigated in this analysis, together with the corresponding S^^^^ . In 
the last column, we give the reference to recent works that employed the given weight function in 
analyses of r decay data. The kinematical weight function w-r was used many times throughout the 
literature and we refrain from quoting all the works that employed it. 

in the existent analyses was mainly the enhancement or the suppression of condensates and 
DV contributions. In combined fits to sets of moments (e.g. refs. [5-9]), the final value of 
receives contributions from the perturbative terms of all moments employed. Therefore, in 
order to achieve a trustworthy determination of the coupling, it is desirable to understand 
the convergence properties of the perturbative component for all the moments employed in 
the Us analysis. 

In the remainder of this section we study the behaviour of the term b^^(m^^ for the 
moments of table 2 in FOPT and CIPT, given by eqs. (2.13) and (2.15) respectively. In 
doing so, we employ for the coefficient cs^i the estimate of eq. (2.11). Regarding higher 
orders, two scenarios are considered: the reference model of ref. [18], given in eq. (3.6), and 
the alternative model of eq. (3.9), which provides an example case where CIPT is better 
than FOPT for the kinematical weight. In plots for the perturbative series we display the 
results for both models side by side to facilitate the comparison. The respective Borel 
resummed values are also shown, together with the Borel ambiguity. Since both models 
are matched to the first five coefficients, their results for CIPT and FOPT are identical by 
construction up to the fifth order. 

In our analysis, it becomes clear that one can group the 17 different into four classes: 
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the monomial terms, the pinched-weights with a "1" in the weight function, pinched weights 
without a "1", and all moments that contain the term x, which form a separate category. 
We analyse these classes in the remainder of this section. 

4.1 Building blocks: the monomial terms 

Since the weight functions employed are all polynomial, it is instrumental to start the 
analysis with the monomial terms given in the first 5 entries of table 2. 

The first of them, the constant = 1, is quite particular. On the experimental side 
this moment gives the same weight to the whole spectrum. On the theory side, it obtains 
very little contribution from condensates, since no powers of x appear, picking up only the 
logarithmic contributions in the Wilson coefficients of the OPE.^ The behaviour of 6w} for 
FOPT and CIPT is shown in figure 1(a) for the reference model, and in figure 1(b) for the 
alternative model. A feature that can be observed in FO is a rapid growth in the first few 
terms, followed by a decrease in the value of 6w}. In lower orders, the series overshoots the 
Borel summed values for both models. Later, FOPT oscillates around the Borel sum. The 
amplitude of the oscillations are much smaller in the RM, and from order n = 5 the series 
can be considered a good approximation to the true result in the sense of an asymptotic 
series. The smaller oscillations in the RM model are due to the previously mentioned 
cancellations among the Cn and the g'rT^^ series. CIPT is more stable in both cases, but it 
fails to give a good approximation to the Borel resummed value of the RM. In the case of 
the alternative model, CIPT is able to approach the true value, as expected. The smallness 
of the ambiguity of the Borel integral due to the poles on the integration contour (indicated 
by the horizontal shaded band in the figures) in these cases can be understood since the 
moment receives only small logarithmic contributions from the condensates in the OPE (for 
lack of powers of x). Accordingly, the ambiguities of the IR poles are also small. ^ 

The behaviour of 6^} for the monomials = x'^ , W4 = x^ , = x^ is qualitatively very 
similar. Therefore, we display only the representative case of in figures 1(e) and 1(f) for 
the RM and AM, respectively, which highlights these similarities. We note that the values 
of 5w} are 4 to 6 times smaller than the ones for wi. This plays an important role in the 
case of moments with pinching. Finally, ws is maximally sensitive to the Borel ambiguity 
of the IR singularity at u = 3, present in both models. However, the residue in the case of 
the RM is about 5 times smaller than in the AM, which explains the different magnitudes 
of the shaded bands in the two plots. 

The behaviour of the monomial ^2 = 2; is exceptional, as shown in figure 1(c) for 
the reference model, and in figure 1(d) for the alternative model. One can separate the 
behaviour in two parts. In the first terms, the FOPT series again grows rapidly and then 
decreases. This is a common feature in the other monomials as well. For higher orders, in 
the RM, FOPT never reaches a plateau: Sw^ exhibits "run-away" behaviour and decreases 
monotonically from the 3rd order. The FOPT series shows no sign of stabilisation around 

^This is precisely the reason why this moment is central to the analysis of refs. [8, 9], where one wants 
to extract the DVs from data. 

^They would be zero if the poles were simple. Since we include the four-loop structure for the renormalon 
singularities, they are small, but non-zero. See Appendix A of ref. [18] for the explicit formulae. 
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Figure 1. Sw} for wi — 1, W2 — x, and W3 — x^, as a function of the order up to which the 
perturbative series are summed for FOPT (black) and CIPT (gray). The horizontal bands give the 
Borel resummed result. The left-hand figures are for the RM of ref. [18], the right-hand ones are 
for the alternative model of eq. (3.9). We use as{mr) = 0.3186. 



the true value, though it develops an inflection point close to the Borel sum. The sign of 
the run-away behaviour (negative) is correlated with the sign of the x monomial (positive) . 
For the alternative model, FOPT still oscillates around the true value. CIPT, on the other 
hand, is rather stable until the onset of asymptoticity, and provides a good approximation 
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in both models. The results in this prescription are within the Borel resummed values for 
the RM, given the larger ambiguity in this case. The larger ambiguity stems from the fact 
that the moment is maximally sensitive to the gluon condensate contribution. This gives a 
larger contribution from the ambiguity of the IR singularity at ti = 2, which is (artificially) 
absent in the AM. 

From this discussion of monomial moments we can already extract a few important 
observations: 

• In the AM CIPT always provides the better approximation. This seems to be a 
generic feature of the AM, not restricted to the kinematical weight. 

• In the reference model, the monomial ^2 = is problematic for FOPT due to run- 
away behaviour, which is correlated with the large d = 4 condensate contribution to 
this moment. 

• For the other monomials low-order approximations in both FOPT and CIPT are 
problematic in the RM case. However, while FOPT converges to the Borel sum for 
n ^ 5, CIPT never reaches it. This is similar to the behaviour found in ref. [18] for 
the kinematical weight. 

In the following, we discuss the remaining moments wq to wn which are composed of the 
monomial terms. Their behaviour can essentially be understood as a linear combination of 
what has been discussed in this section. 

4.2 Pinched weights with a "1" 

We now turn our attention to moments with pinching and start with moments that contain 
a term "1" in the weight function and that do not have a linear term x. In table 2, these 
moments are wj, Ws, Wr- For all these moments, the term 1 sets the scale and the higher 
powers only introduce corrections to this leading result. Since they are pinched moments, 
they always have at least one negative term. This leads, in general, to a stabilisation of 
FOPT in both models. 

In the case of FOPT for the RM, the monomials in Wr and u^y conspire to give an 
excellent cancellation of the initial overshooting, leading to a series that approaches the 
Borel sum very fast. The cancellations between CnS^^'^ and gl^'"'^ then make FOPT rather 
stable after approaching the true value. This can be observed in the results for the RM 
given in figures 2(a), 2 (c) , and 2 (e) . The results for tfg = 1 — resemble more the ones for 
wi = 1 since the corrections due to are quite small. CIPT, on the other hand, misses 
the cancellations and never approaches the Borel resummed values. We also observe that 
the CIPT series enters the sign alternating regime earlier than FOPT. The corresponding 
results for the AM are shown in figures 2(b), 2(d), and 2(f). As foreseen, here CIPT tends 
to give a better approximation to the resummed series. Albeit more stable than in the case 
of the monomials, FOPT still displays oscillations around the Borel sum. 

There are many other possible moments, not shown in table 2, that display a very 
similar behaviour. We have investigated the family w{x]a) = 1 + ax"^ — {a + l)x^ for several 
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Figure 2. 5^' for Wr, w^, and wg, as a function of the order up to which the perturbative series 
are summed for FOPT (black) and CIPT (gray). The horizontal bands give the Borel resummed 
result. The left-hand figures are for the RM of ref. [18], the right-hand ones are for the alternative 
model of eq. (3.9). We use as{mr) = 0.3186. 



different values of a witli results rather similar to the ones discussed above. Other moments 
that do not start with the unity, but with another constant of the same order such as 
w{x) = |(1 — x)^(l -|- x)(l -|- a; -|- Ax"^) also give qualitatively similar results. Thus, the main 
observation for this class is that 
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• pinched-moment weights with a "1" but no hnear term all behave similar to the kine- 
matic weight Wr, favouring FOPT over CIPT for the reference model and vice versa 
for the alternative model. In each case the perturbative expansion in low orders 
approaches the Borel sum rather quickly. 

4.3 Pinched weights without a "1" 

The next class of moments that we analyse here are moments with pinching but that do 
neither contain a constant term nor one linear in x. As examples for this group we employ 
Wi6 and Wi7 of table 4; two of the triply-pinched moments used by ALEPH and OPAL [3-6]. 
In ALEPH's notation, these moments are part of a family denoted by w^^'^\ and read 

^^(i.'^) = (1 - xfx''{l + 2x) = x^ - x^^^ - 3 + 5 - 2 x^^^ . (4.2) 

In table 4 we have w\-i = w^^'^'^ ^ = w^^'^\ wiq = w^^''^\ and wn = w^^'^\ Besides 
the "1" in wi^, both the first two also contain a linear term in x and will be discussed in 
the next section. From table 2 we see that the pinched weights without a "1" have very 
small 6^^^ but 0(1) coefficients of the monomials. This enhances the relative importance 
of power corrections. 

We consider first the case of wiq in the RM, shown in figure 3(a). As expected, 6wie 
tiny, almost 50 times smaller than the corresponding correction for wi = 1. The combination 
of powers of x does not improve the behaviour of FOPT that overshoots largely the Borel 
sum in the first few orders. It eventually approaches the Borel result for higher orders, 
just before the onset of asymptoticity. The bad behaviour of CIPT already observed for 
the monomials is amplified and the CIPT series goes astray. The situation in the AM is 
somewhat improved, but mainly due to the large Borel ambiguity associated with the IR 
pole at It = 2, see figure 3(b). Again, FOPT displays large oscillations around the Borel 
result, whereas CIPT grows monotonically away from the resummed result before the sign- 
alternating asymptotic behaviour sets in. As seen in figure 3(c), for FOPT in the RM, the 
perturbative contribution to the moment of wn has a behaviour qualitatively similar to 
wiQ, though for higher orders it is slightly more stable. Also CIPT approaches the Borel 
sum before the series becomes asymptotic after the 9th order. In the AM, figure 3(d), 
both, FOPT and CIPT fail to approach the Borel sum, FOPT once more displaying large 
oscillations.^ We therefore conclude: 

• The perturbative expansions for this class of moments tends to be unreliable in both 
FOPT and CIPT, and independent of the model for the unknown higher-order coef- 
ficients. 

* A possible criticism against our analysis of wn within the reference model could regard the lack of an 
IR singularity at u = 4. Since the moment starts with it is maximally sensitive to D = 8 contributions 
in the OPE, which corresponds to the ambiguity of the IR renormalon at u = 4. We investigated this issue 
by considering a model where one adds an IR renormalon at it = 4 and leaves only a constant dp ° in the 
model. (This model is briefly discussed on page 24 of ref. [18].) After performing the matching, the residue 
of the renormalon at m = 4 turns out to be small (d^^ = 5.64) and the changes in the other parameters 
negligible. The additional Borel ambiguity arising from it = 4 is also small. Therefore, the result shown in 
figure 3(c) is not altered in any significant way, which corroborates the assumption that the singularities at 
u = 2 and u — 3 are the dominant ones. 
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Figure 3. 5^} for wig and wn, as a function of the order up to which the perturbative series 
are summed for FOPT (black) and CIPT (gray) are summed. The horizontal bands give the Borel 
resummed result. The left-hand figures are for the RM of ref. [18], the right-hand ones are for the 
alternative model of eq. (3.9). We use asimr) = 0.3186. 



• Pinched moments without a "1" are sensitive to condensates, but the poor perturbative 
approximations render condensate determinations from these moments unrehable. 
This conclusion appears to be largely model-independent. 

4.4 Moments containing a term x 

We have relegated to this section the analysis of weight functions containing the mono- 
mial X. This choice is based on the observation that the behaviour of these moments, which 
are maximally sensitive to the D = A correction in the OPE, is qualitatively different in 
the RM of ref. [18]. This was shown for the monomial above and here we discuss pinched 
weights containing this term. There are several of them in table 4: tfg, wg, wiq^ wh, ifis, 
i(;i4, and wi^. Again, they display very similar qualitative behaviours and it suffices to 
expose in detail only three representative examples. 

We start with the simple case oi = 1 — figure 4(a). The term 1 sets the scale, but 
now, since the perturbative series for the monomial W2 = x decreases monotonically, the 
perturbative series for wg, whose linear coefficient has negative sign, grows monotonically. 
The result for FOPT crosses the Borel resummed value around the 7th order where it also 
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Figure 4. S^. for wg, W13, and W15, as a function of the order up to which the perturbative series 
are summed for FOPT (black) and CIPT (gray). The horizontal bands give the Borel resummed 
result. The left-hand figures are for the RM of ref. [18], the right-hand ones are for the alternative 
model of eq. (3.9). We use as{mr) = 0.3186. 



has an inflection point. CIPT can only approach the Borel sum shortly before it becomes 
asymptotic around the 9th order. The situation is very similar when higher orders of x are 
added to the weight function. In figure 4(c), the moment 1^13 is displayed as a representative 
example. (It corresponds to w^^'^^ in ALEPH's notation [3-5] of eq. (4.2)). This weight 
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starts with 1 — x followed by higher-order terms in x. Qualitatively, the only difference now 
is that CIPT never comes close to the Borel sum. The picture changes in the case of W15 
(or t/;^^'^^) since the term 1 is missing. As figure 4(e) shows, this case is more similar to 
the monomial x itself, though the higher powers in x soften the behaviour. For CIPT, the 
situation is not much different either, but it only approximately represents the Borel sum 
at low orders. 

For the alternative model which does not contain the renormalon singularity at u = 2, 
the first two moments we and wis, figures 4(b) and 4(d), are similar to what was observed 
in the case of pinched moments with the term 1. Now, however, CIPT approaches the Borel 
results less fast in the case of wq, and not at all in the case of wis- FOPT in both cases 
displays oscillations around the true value. The last case, that of wi^ shown in figure 4(f), 
is less satisfactory. As was the case for other moments from weight functions starting with 
a power of x, the values of 6w^^ are small. While CIPT misses the Borel sum completely, 
FOPT only approaches it around the 10th order. We summarise our observations as: 

• Weights without a "1" are again unreliable in FOPT and CIPT, and in both models, 
especially at intermediate orders. 

• In the RM, FOPT exhibits run-away behaviour and CIPT may not approach the 
resummed result. Overall, the perturbative expansion does not behave as well as 
for moments without a linear term, which is related to the sizeable = 4 power 
correction. 

• In the AM there is no clear preference for one of the two methods. 
4.5 Main lessons from the moment analysis 

While we already summarised our main observations for each class of weight functions, we 
collect again here the most important points. 

Some of the pinched-moments (with the "1", without the "x") display a particularly 
fast convergence of FOPT towards the Borel resummed values, especially the moments Wr 
and wj, as shown in figures 2(a) and 2(c). On the contrary, in Borel models that contain a 
u = 2 pole residue of natural size, CIPT generally does not approach the Borel sum before 
the divergence of the series sets in. This is different in the alternative model, where the 
u = 2 pole is artificially suppressed, and it coincides with the main findings of ref. [18]. 
Thus, if the reference model is adopted as the most plausible one (as we would do), one 
again arrives at the conclusion that FOPT provides a better approximation than CIPT, 
also at order n = 4, 5. 

The investigation of moments also reveals that the qualitative behaviour of their per- 
turbative expansion depends only on a few features of the moment and the model of the 
Adler function. As concerns the model, we have already emphasised the crucial question 
of the size of the residue of the u = 2 singularity that corresponds to the D = 4 power 
correction, which motivated the choice of the AM. 

As concerns the moment function itself, an important observation is that moments 
that start with high powers of x, such as W15 and wiq, employed by the ALEPH and OPAL 
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collaborations, have a bad behaviour of the perturbative series. For these moments, neither 
FOPT nor CIPT are able to provide a decent approximation to the Borel resummed values 
in the first few orders. Also, these moments have a very small value of 6^} , with large 
relative Borel ambiguities, which makes the reliable separation of power corrections from 
the uncertain perturbative approximation problematic. Therefore, these moments are not 
an optimal choice for an as analysis. (This has already been pointed out in ref. [7].) 

Finally, for moments that contain a linear term x, the reference model, or others that 
include an IR pole at m = 2, both FOPT and CIPT behave badly. In the case of FOPT, 
the series is quite unstable, which results in large errors due to the truncation of the series, 
producing unstable results for Og. (This was noticed — in practice — in the exploratory 
fits of ref. [40], and it is the reason why, in refs. [8, 9], moments with the term x were not 
considered.) The situation of CIPT for these moments is also unsatisfactory because the 
series are unstable and/or do not approach the Borel resummed result. This suggests that 
moments with a linear term x should also be avoided in Og determinations. 

5 Validation of the reference model 

In the previous section we learnt that for moments with good perturbative convergence the 
comparison of FOPT and CIPT leads to the same conclusion as for the inclusive hadronic 
tau width studied in ref. [18]. Hence, the crucial factor in deciding whether FOPT or CIPT 
should be the method of choice remains the plausibility of the reference ansatz for the Adler 
function (favouring FOPT) as compared to, e.g., the alternative model (favouring CIPT). In 
addition to the general arguments for the reference model reviewed in section 3, we discuss 
in this section two further checks, one inspired by ref. [24], which support the plausibility 
of the ansatz and results of ref. [18]. 

5.1 Adding a polynomial term 

In ref. [18], the known higher-order coefficients, plus an estimate for cs^i, are used to fix the 
residua of the renormalon poles, while the first two polynomial terms are obtained by also 
fitting the coefficients ci^i and C2,i. It is assumed, therefore, that the renormalons dominate 
at intermediate (and higher) perturbative orders. However, the fact that is small in 
the reference model of [18], indicates that C2,i is already well saturated by the renormalon 
poles. The procedure has been criticised in ref. [24], where the authors argue that the 
truncation of the polynomial terms at linear order is arbitrary. They propose to add a 
term to the polynomial and study the behaviour of those models when the coefficient 
is fixed to six different values: df^ = -1, -0.5, 0, 0.25, 0.5, 1. For the value = 0, the 
RM is recovered. 

The inclusion of a fixed term in the modelling of the Borel transform of the Adler 
function has consequences for the residua of the renormalon poles. They have to adjust to 
the existence of this term which contributes to cs^i, which leads to abnormally high values 
for the residue of the IR pole at n = 3, see table 3. Consequently, large cancellations among 
the contributions of the IR poles sX, u = 2 and u = 3 arise. In the two extreme cases studied 
in [24] , (ig*^ = il, the residue of the pole at « = 3 changes by factors of —11 and 13 with 
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dr 


dPO 


df o (fixed) 


RM [18] 


3.16 


-13.5 


-1.56 X 10-2 


0.78 


7.66 X 10-3 





[24] 


-1.50 


149.7 


-5.90 X 10-2 


10.68 


3.85 


1 


-u^ [24] 


7.82 


-176.8 


2.79 X 10-2 


-9.12 


-3.83 


-1 



Table 3. Parameter values of three models for the physical Adler function. "RM" represents the 
central reference model of ref. [18]. The models denoted by "±u^" are those discussed in ref. [24] 
where the polynomial coefficient ^2*"^ is taken to be ±1. 



respect to the RM of [18], for = +1 and = —1 respectively. Furthermore, fixing 
the v?' term also forces a break-down of the renormalon dominance of the coefficients C2,i 
and cs^i. This is apparent from the values of the other polynomial terms. In the extreme 
cases d^^ = ±1, the module of the coefficient dg*-* is more than 10 times larger than in 
the RM. The next coefficient, df*-*, is also much larger, more than 100 times the one found 
in [18]. The values for the residua and the polynomial terms for the RM and for the extreme 
cases d2^ = ±1 are given in table 3. Finally, the coefficient d^'^ and the u = 2 residue d2^ 
share an almost linear relation, such that d^ vanishes for d^^ = 0.678. This particular 
case constitutes another model for which CIPT generally better approximates the Borel 
sum for 5w}- 

The inspection of the Adler function D{s) on the complex circle s = m^e^'^ sheds further 
light on the plausibility of the models considered here. It is expected that the perturbative 
expansion breaks down in the vicinity of the physical, Minkowskian axis ((/> ~ or (/) ~ 27r), 
but that it should work well in the Euclidean region (f) it. The behaviour of Yie[D{(l))] 
along the complex contour is displayed in figure 5, of which the upper plot corresponds 
to FOPT and the lower to CIPT. The dotted, dot-dashed and double-dot-dashed curves 
are the 3rd, 4th and 5th order purely perturbative results respectively. The thick solid 
line corresponds to the Borel sum of the reference model. An analogous plot was already 
shown as figure 9 in Appendix B of [18]. In addition, we now also display the additional 
models with d2^ = 1 (short-dashed line) and d2*^ = — 1 (long-dashed line). Furthermore, 
the shaded area indicates the 5th order PT result when the coefficient cs^i is varied in the 
range 05,1 = 283 ± 283. 

The following observations can be made on the basis of figure 5. For an asymptotic 
expansion, the last included term should provide an approximate error estimate for the full 
sum. Though this is strictly true only for sign- alternating asymptotic series, we expect the 
estimate not to be wildly violated. Employing the shaded area as such an error estimate, 
it is seen that in the Euclidean, ~ vr, the RM of [18] lies rather close to this region. 
On the other hand, the models with d2^ = il, even in the Euclidean domain where PT 
should work well, lie far from 5th order perturbation theory. Furthermore, moving away 
from the Euclidean axis, strong oscillations in Re[Z)((/>)] are found in those models. It seems 
rather unlikely to us that QCD behaves in this way. Turning the argument around and 
investigating which values of d2 ^ would yield models compatible with the shaded area, we 
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Figure 5. Real part of the Adler function b{s) on the complex circle s = m^e"^. Borel sums 
are displayed for the RM (solid line), ^2'"^ = 1 (short-dashed line) and 1^2'^ = — 1 (long-dashed 
line). The dotted, dot-dashed and double-dot-dashed lines correspond to 3rd, 4th and 5th order 
of perturbation theory. Finally, the shaded area indicates the 5th order PT result while varying 
C5,i = 283 ± 283. The upper plot shows FOPT and the lower CIPT. We employ as{mr) = 0.3186. 



roughly obtain the range —0.55 < ^2'^ < 0. 

To summarise, there are two arguments in favour of the procedure adopted in ref. [18] 
for the treatment of the polynomial terms in the ansatz of eq. (3.6). First, the fact that 
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turns out to be so small in the central model fit, together with the observed hierarchy 
cIq^ ^ di^j which leads to a renormalon dominance of the coefficients at orders as low as 
a^. Second, the unnaturalness of the Adler function shape along the circle when values 
of 1^2^! are imposed in the extended model suggested in ref. [24] that seems to obstruct 
duality even in the Euclidean region. 

5.2 Matching in the large-/3o hmit 

Another check whether a simple ansatz such as eq. (3.6) can work can be derived from the 
large- /3o limit. As discussed in section 3, an analytic result for the Borel-transformed Adler 
function is available in this limit, eq. (3.4), and hence the exact perturbative coefficients 
Cn,i are known to all orders [36, 37]. Here, we propose to emulate the matching procedure 
performed in QCD in the context of the large-/3o approximation. That is, we make a simple 
ansatz similar to the RM and fit the parameters of this ansatz to the low-order Cn,i in the 
large-/3o approximation. We then compare the so-obtained model for the higher-order terms 
to the exactly known ones. 

In order to implement the matching procedure in the case of the large- /3o limit, we 
adapt the reference model to the present case by using simple and double poles (instead of 
branch cuts) for the renormalon singularities. In the spirit of eq. (3.6), the new model can 
then be written as 

^l^K") = ^ + + (if^ + 4° + ° " + (5., 

When emulating the procedure of ref. [18], we set the term = 0. In large-/3o, the IR 
pole at u = 3 is a double pole and therefore 73 = 1. We simulate our ignorance of the 
structure of this pole in full QCD by taking either 73 = or 73 = 1. From eq. (5.1), by 
changing 73 and the assumptions about (ig*^, we define five different models and perform 
the matching to the first coefficients of the exact large-/3o Adler function. The models will 
be denoted by the values of these parameters as M(73;d2^). In the remainder of this 
section the characteristics of these models are discussed, the matching is described, and we 
compare with the exact large- /3o limit. 

We begin by performing the matching treating the term ^2^ as suggested in ref. [18], 
which is equivalent to setting ^2^ = 0. Furthermore we employ 73 = 1, which gives a 
double pole for the IR singularity at n = 3, in agreement with the exact result eq. (3.4). 
This model is referred to as M(1;0). The renormalon residua are fixed to the coefficients 
C3,i; C4^i, and cs^i; the two polynomial terms are found by enforcing the true large-/?o values 
of ci^i and C2,i. The results of this model are shown in the second row of table 4. As 
in the full QCD case, the value of df*-* turns out to be small, and significantly smaller 
then cIq*"^, indicating the renormalon dominance at intermediate orders. Furthermore, the 
residua of the renormalon poles are in the ball-park of the true results, although they have 
to compensate for the lack of higher order poles. 

Since the purpose of the model introduced in ref. [18] was to decide upon the best way 
to perform the RG improvement of the series, it is legitimate to ask if the description of 
higher orders is successful in the present case. This can now be unambiguously tested by 
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1 




U 


1 


4° 


1 PIT CTf^— /"In 


1 7 S4 


— 1 n 4Q 




-2 








M(1;0) 


16.53 


-45.79 


4.10 X 10" 


-2 


-2.90 


-0.44 


(fixed) 


M(0,0) 


8.34 


-18.43 


4.46 X 10" 


-2 


2.25 


0.27 


(fixed) 


M(l;free) 


18.85 


-55.02 


3.92 X 10" 


-2 


-3.03 


-0.34 


5.79 X 10-2 


M(l,l) 


56.54 


-205.27 


9.76 X 10- 


-3 


-5.15 


1.31 


1 (fixed) 


M(l,-1) 


-23.46 


113.67 


7.23 X 10" 


-2 


-6.52 X 10"^ 


-2.19 


-1 (fixed) 



Table 4. Residues of poles and polynomial parameters of the models discussed in the text. The 
first row gives the exact result in the large-/3o approximation. The models are defined by the values 
of 73 and ^2 *^ in eq. (5.1) and denoted by M{-f3; d^'^). 





C6,l 


C7,l 


C8,l 


C9,l 


ClO.l 


Cll.l 


large-^o 


-1.99 X 10^ 


9.86 X 10'' 


-1.08 X 10^ 


2.78 X 10^ 


-5.39 X 10^ 


1.40 X IQi" 


M(1;0) 


-2.46 X 10^ 


1.08 X 10^ 


-1.30 X 10*^ 


3.28 X 10^ 


-6.68 X 10^ 


1.74 X IQi" 


M(0;0) 


-3.53 X 10^ 


1.08 X lO-"^ 


-1.51 X 10^ 


3.46 X 10^ 


-7.38 X 10^ 


1.87 X 10^° 


M(l;free) 


input 


1.07 X 10^ 


-1.21 X 10*^ 


3.17 X 10^ 


-6.35 X 10^ 


1.67 X 10^° 


M(l;l) 


5.69 X 10^ 


8.54 X 10* 


2.74 X 10^ 


1.39 X 10^ 


-9.09 X 10^ 


4.95 X 10^ 


M(l;-1) 


-1.06 X 10^ 


1.30 X 10^ 


-2.88 X 10^ 


5.17 X 10^ 


-1.25 X 10^ 


2.98 X IQi" 



Table 5. Higher-order coefficients from the five models described in the text compared with the 
exact large-/3o results. Models are defined by the values of 73 and d^'^ in eq. (5.1) as M(73;(i2°). 

comparing the higher order coefficients predicted by the fitted Adler function with the exact 
results in large-^So- Numerically, this comparison is shown in table 5. Graphically, FOPT, 
CIPT, and the resummed results for the model are sliown in figure 6(b) (for Wt), whereas 
the equivalent plot for the exact large-/3o results is given in figure 6(a). The qualitative 
agreement of the results of table 5, together with the striking similarities of figures 6(a) 
and 6(b), demonstrate that — in spite of the simplifications of the model with respect to 
the exact results — the model reproduces faithfully the FOPT and CIPT series up to higher 
orders, as well as the Borel resummed value. 

In the previous model, the parameter 73 was fixed to the true value of the large- /3o limit. 
Let us investigate the consequences of a wrong choice for 73. For full QCD this is a relevant 
open issue as several operators contribute at D = 6, and the general renormalon structure 
has not yet been established. We define a new model that differs from the previous one 
only by having a simple pole at u = 3, thus 73 = 0. Accordingly, we denote this model by 
M(0; 0). Performing the matching, we find the results of the 3rd row of table 4. The choice 
73 = enforces a number of adjustments in the parameters of the model with respect to 
the case where 73 = 1. The residue of the pole at u = 2 changes and is less well reproduced 
than in the previous case. The polynomial terms are different, but the hierarchy between 
the terms is still preserved. In particular, inspection of the third row of table 5, together 
with figure 6(c) shows that the higher-order behaviour of FOPT and CIPT continues to 



-25- 



0.28 
0.26 
0.24 
0.22 
0.20 
0.18 
0.16 
0.14 
0.12 



0.28 
0.26 
0.24 
0.22 
0.20 
0.18 
0.16 
0.14 
0.12 



T ' 1 ' 1 ' r 

FOPT — • — 
CIPI 
Borel sum 




(a) Exact large-/?o 



11 



13 



15 



1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 
FOPT — •— 
CIPT 

Borel sum , 




• 








/ \ 










• / 

1 / . 1 . 1 . 1 , 1 , 


1 , 1 





15 



0.28 
0.26 
0.24 
0.22 
0.20 
0.18 
0.16 
0.14 
0.12 



0.28 
0.26 
0.24 
0.22 
0.20 
0.18 
0.16 
0.14 
0.12 



1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 
FOPT — •— 
CIPT 

Borel sum %. 

~ * / \ 






















• / 

1 / , 1 . 1 , 1 , 1 , 


'.I , 1 





11 



(b) M(1;0) 
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Figure 6. Values of Swliiri^) as a function of the order n up to which the perturbative series has 
been summed for FOPT (black) and CIPT (gray). Horizontal gray lines give the Borel resummed 
result and the bands give the estimated ambiguity. In (a) one sees the result from the exact large-/3o 
limit. In (b)-(f), we show results from models matched to the first few coefficients of the large-/3o 
Adler function (see text). Models are defined by the values of 73 and in eq. (5.1) as M(73; ^2^^). 
For consistency with large-/3o, we perform the ctg running at one loop. We use as{m'^) = 0.3186. 
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resemble very much the exact result. In agreement with figure 6(a), FOPT approaches 
the Borel resummed value better than CIPT. Still, one must admit that there is a shift 
in the Borel resummed value of ~ 5%, which is beyond the ambiguity of the exact result 
(~ 1%). Nevertheless, the conclusions about FOPT versus CIPT — and the superiority of 
the former — remain intact. 

We now investigate the question of adding a term d^'^u'^ in eq. (5.1). Since all co- 
efficients Cn,i of the Adler function are available in the large-/3o approximation, we have 
the freedom to include the coefficient cg,! in the matching. This allows us to keep the six 
parameters of the model free, including ^2^, and then follow a strategy similar to the one 
above (we use 73 = 1). This model is termed M(l;free). We employ the coefficients 04,1, 
C5,i, and C6,i to fix the residua of the renormalon poles, while the polynomial terms are fixed 
by C3^i, C2,i, and ci^i. The results after matching are found in the third row of table 4. The 
inclusion of the free parameter preserves the hierarchy ^ df*-* S> (^2^, indicating 
the renormalon dominance at intermediate and higher orders. This is in agreement with 
what was found before, where we kept d^'^ = fixed. Table 5 and Fig 6(d) show that the 
quality of the description of higher orders remains impressive when we introduce the term 
^2*^ keeping it free and using the procedure of ref. [18] to fix it. 

Let us turn now to the consequences of having a fixed parameter d^'^ . To make contact 
with ref. [24], we study the two cases ^2^ it 1 (again 73 = 1). These configurations are 
denoted M(l,±l). In table 4, one sees some issues that emanate from this choice. The 
hierarchy of the polynomial terms is broken and, what is more, the values of the residua 
of the IR renormalon poles are very different from the exact ones. We observe large values 
of the residua, leading to strong cancellations between the contributions of the different 
IR poles, a feature that was already observed in full QCD. What is then the most faithful 
description of the higher orders? Table 5 demonstrates that the description with fixed 

^ o{l) is very poor. The large-order coefficients are badly reproduced; even the sign 
is wrong in two of them. Figures 6(e) and 6(f) show that the perturbative series goes astray, 
and the Borel resummed results are very different from the exact one shown in figure 6(a). 

This exercise shows that the model is incompatible with a v? term whose coefficient is 
of order unity. In real QCD, only four coefficients of the Adler function are known exactly. 
In such a scenario, it seems that the best strategy is to keep ^2^ = 0, since we learn in 
large-/3o that fixing this parameter to an arbitrary value is not a better option. A final 
comment is order. Figure 6 shows results for the kinematic moment Wr only. We have 
studied all moments displayed in table 2 also for large-/3o and the conclusions regarding the 
quality of the description are not altered in other cases. 

6 Consequences for the determination of cxs 

Based on the perturbative behaviour of the Adler function under two different assump- 
tions for the higher-order coefficients, we argued in section 4 that some weight functions 
are more suitable for an as analysis from hadronic r decays than others. The aim of this 
section is to corroborate these findings comparing the predicted moments with experimen- 
tal values. We want to check the internal consistency of the predictions from different 
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moments and study the uncertainties associated with them once all the contributions are 
taken into account (power corrections and duality violations). Rather than comparing mo- 
ment predictions with data for some given value of a^, we perform this consistency check 
by determining the values of Os from each individual moment. The spread of values 
and uncertainties then provides the desired information. We emphasize that the aim of this 
exercise is not a precise determination of the strong coupling. 

We consider FESRs for the weight-functions of table 2. In the notation introduced in 
eq. (2.6), we study the moments i?y'^^(m^), relying on external input for the power correc- 
tions and for DVs. A more comprehensive analysis should include more than one moment 
and treat consistently all the ingredients of the theoretical description: the perturbative 
series, power corrections, and duality violations. In such an analysis, an involved non-linear 
multi-parameter fit taking into account all correlations is unavoidable (see e.g. [8, 9]). Here 
our error bars in as are smaller than the ones of a self-consistent analysis due to the fact 
that we do not perform a multi-parameter fit. 

The different theoretical components in the computation of the moments are treated 
as follows. As discussed in section 2, the perturbative FO and CI series are summed up to 
order using cs^i = 283. The uncertainty due to the truncation of the series is estimated 
by either taking cs^i = or C5.1 = 566. A third as value is obtained from the Borel 
resummed results of the reference model. 

The power corrections in the OPE are taken as an external input. Corrections due to 
dimension 4, 6, and 8 are considered. (The mass corrections, D = 2, are also taken into 
account although they can safely be neglected due to the smallness of the quark masses.) 
For = 4, we use the gluon condensate value (aG^) = (0.012 it 0.012) GeV^, and include in 
the coefficient function the known as corrections with their logarithms [41]. At dimenstion 
six, the main contributions arise from the four-quark condensates, since the coefficient of the 
three-gluon condensate vanishes at leading order and all terms proportional to quark masses 
can safely be neglected. One usually resorts to the vacuum saturation approximation to 
write these contributions in terms of squares of the quark condensate [44], introducing the 
parameter pv+A to account for deviations from this assumption. The D = 6 corrections 
are then proportional to pv+a{qq)'^ ■ In our estimates we use pv+A = 2 it 1 [18] and 
{qq){mr) = —(272 ± 15 MeV)^ [45]. A crude estimate of = 8 is included adding the a 
term Cgy+A/ s'^ to the Adler function. To evaluate the impact of this contribution we use 
Cgy+A = (0 ± 5) • 10^'^. Our estimates for D = Q and 8 agree, within uncertainties, with 
the results of the fits of ref. [9]. 

A phenomenological estimate of the longitudinal contribution from scalar and pseu- 
doscalar correlators is included in the spirit of refs. [42, 43]. The weight-function depen- 
dence of the non-perturbative corrections make some of the moments quite insensitive to 
the details of the non-perturbative contributions. However, for other weight-functions, this 
is not the case and the final as values depend heavily on the non-perturbative input. 

Finally, the contribution from DVs to the moments is computed using the results of 
ref. [9]. The corresponding term takes the form 




(6.1) 
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where the DV part of the spectral functions, is parametrised as [15] 

Pv/Ais) = exp {-8v/A - 1v/as) sin {ay /a + I^v/as) ■ (6.2) 

The eight parameters for the description of DVs are taken from the results of a fit to y 
and A using updated OPAL data shown in table V of ref. [9] (we take the FOPT fit with 
■Smin = 1-5 GeV^). Their correlations are employed in the Monte Carlo estimate of the error 
induced by DVs. Note, however, that within our simplified Ug determination in which DVs 
and power corrections are not determined self-consistently, the precise DV parameters are 
not important for the following. Within other errors, the DV term (6.1) has a negligible 
impact on from the moments discussed below. 

The experimental counterparts of the moments i?y'^^(m^) are obtained performing a 
discretised version of the integral given in eq. (2.2), where the spectral functions are the 
ones from the ALEPH collaboration after the 2008 update performed in ref. [5]. It is known 
that in this data set contributions to the correlations due to the unfolding procedure were 
inadvertently omitted [40]. At present, in a precise determination of a^, the safest option 
is to turn to an updated version of the OPAL data [6, 9]. Nevertheless, for our exploratory 
purposes, this omission is harmless and would only alter the experimental error bars, that 
are potentially underestimated. 

In figure 7, we display results of this simplified analysis for 13 of the weight func- 
tions of table 2, and using FOPT, CIPT, and the Borel sum within the RM. (We do not 
show results for the monomials W2 to w^.) The inner error bars give the experimental er- 
rors, while the external error bars include the theoretical error as well. The relative size of 
the experimental uncertainties depends strongly on the moment considered. This is under- 
standable since the spectral functions have larger relative uncertainties in the higher-energy 
part. Moments that emphasise this region (or that do not suppress it) such as wi{x) = 1 are 
penalised and have significantly larger uncertainties. In general, pinched moments suppress 
the edge of the spectrum and have smaller relative uncertainties. 

We consider several sources of theoretical errors in a^. The first one is the truncation 
of perturbation theory, which is estimated by varying the coefficient Cs^i = 283 it 283. 
Another source of theoretical error is the residual renormalisation-scale dependence. To 
estimate this uncertainty, we re-express the FO series in terms of a different scale rather 
than using a(m^), and vary this scale. The residual dependence is estimated in CIPT 
in a similar fashion by setting the scale to /i^ = —^"^sqx in eq. (2.12). This generates 
additional logarithms that must be taken into account in S^S) . In the Borel resummed 
model, this scale/model uncertainty is also estimated by taking = ^^m^, determining 
a(^mT-), and evolving the result back to ^ = m,-. In all cases, the scale is varied in the 
interval 0.5 < — 1-5 m^. The uncertainty due to D = 4 contributions is estimated 
varying the gluon condensate in the interval (oG^) = (0.012 it 0.012) GeV^. Uncertainties 
from D = Q and D = 8 are computed from the propagation of the error on the quantities 
Pv+A, {qq), and Csy+A- Finally, the uncertainty due to the DV term is estimated from 
a Monte Carlo sample of parameters generated according to the results found in table V 
of ref. [9]. The outer error bars in fig. 7 contain the sum in quadrature of all these errors 
together with the experimental uncertainty. 
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Figure 7. Results of asirrir) from V + A sum rules constructed with the weight functions of table 2 
and using the 2008 version of ALEPH data [5] . The absence of a point means that no reasonable 
value of as was found. Inner error bars give solely the experimental error; outer bars include the 
theory error as well. In the lower shaded band we show explicitly the constant term and the term 
proportional to x for the weight functions (when present). All the moments are pinched except for 
wi. The world average of as is that of the PDG [46]: ^^(m^) = 0.3186 ± 0.0056. 



Let us turn now to an analysis of the results shown in Fig. 7. A first important point is 
that for the moments considered ideal in the study of 6w} , namely, wy, ws, and Wr, we obtain 
an as value compatible with the world average within our simplified analysis. The results 
from FOPT and the reference model are particularly close to the world average, while 
CIPT lies about two sigma away. This is a consequence of the fact that 6w} dominates 
these moments, and the FOPT perturbative series shows a good behaviour. They are 
quite insensitive to the the power corrections and DVs and, accordingly, have relatively 
small theoretical uncertainties. Their experimental uncertainties are smaller than in other 
moments due to the suppression of the edge of the spectral functions. 

The situation is radically different for w^, W15, and wn for which neither FOPT nor 
CIPT were able to give a reasonable value of as- The Borel resummed result from the 
reference model of ref. [18] does not yield any acceptable value for for wiq and wn 
either. For wi4^ and W15 abnormally high values with huge uncertainties are obtained. This 
is a manifestation of the bad convergence properties observed for the family of moments 
defined in eq. (4.2) and used e.g. in refs. ]3-6], combined with the fact that these 
moments receive very large contributions from power corrections (of the order of 50% for 
wiy). Note that with the present treatment of power corrections the theoretical error of these 
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moments is by far dominated by perturbative (scale and c^^i variations) uncertainties. This 
corroborates the conclusion of section 4, namely, that these moments are not the optimal 
choice for an analysis from r decays. For wiq, FOPT and CIPT yield a value of a<j 
with large errors dominated by power corrections of dimension 6 and 8 with a sizeable 
contribution from the perturbative series. The result in FOPT is particularly unstable, 
leading to larger uncertainties. 

The scenario for the moments that contain an x term and the unity lies in between 
the previous two. Moments like wq and also wis still yield reasonable values of Og- On the 
other hand, moments like wg, wiq, wh give lower values of as- This can be understood 
since for these moments the D = 4 correction is increasingly important, due to the higher 
coefficients of the x term. In our estimate, these D = 4 corrections are positive which 
leads to lower values of a^. This is not the case in wi^ due to the D = 6 and D = 8 
contributions arising from the terms — + 5x^ in the weight function, that compensate 
the D = 4 corrections. Since these moments are more sensitive to D = 4, they exhibit 
larger theoretical uncertainties arising predominantly from ^j^.^y/^- 

Finally, from tui = 1, a value of compatible with the world average is also obtained. 
This moment has a much larger experimental uncertainty due to the lack of pinching and 
the corresponding higher contribution from the edge of the spectrum. 

A general observation is that the size of the discrepancy between FO and CI also 
depends strongly on the moment. Moments such as wg, wiq, and u^n tend to minimise 
the discrepancy while the kinematical moment Wr, and to some extent wt = 1 — x'^ and 
Ws = 1 —x^, tend to maximise this difference. We have seen that FOPT is a better approx- 
imation to the Borel resummed results in the reference model; therefore, as values from 
FOPT are systematically closer to the Borel resummed ones than CIPT values. 

7 Conclusions 

In the first part of this work, we have analysed the moment dependent features of the per- 
turbative expansion of the Adler function needed in the theoretical description of hadronic 
r decays. A systematic study of this moment dependence is important since it serves as a 
guide for future analyses of as, and provides further insight on the question whether FOPT 
or CIPT is a better framework for ag extractions. To analyse the higher order behaviour 
of the series we employed two models. The first is the reference model of ref. [18], which 
gives — we believe — a plausible representation of the QCD Adler function. In this model, 
FOPT is the prescription that gives the best approximation to the Borel resummed results. 
In order to assess possible model dependencies in our conclusions, we have also employed 
a model where the IR pole at n = 2, which is the most important one for the perturbative 
behaviour in intermediate orders, is artificially suppressed; this description favours CIPT. 
The behaviour of 6^} for a large collection of moments can then be divided into a small 
number of classes. The general behaviour of the perturbative series can be traced back 
to the singularities of the Borel transformed Adler function and to simple features of the 
weight-functions . 

We have shown that some moments have better perturbative properties than others. 
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Our conclusions are based on the inspection of the perturbative series and on a simphfied 
as extraction from each one of the moments. In particular, polynomial pinched weight- 
functions that contain the unity and do not contain a term proportional to x are found to 
be optimal. This is due to the good behaviour of the perturbative series in these cases, where 
at least one of the methods, FOPT or CIPT, provides a good approximation to the Borel 
resummed results. Another positive feature of these moments is the small contamination 
by power corrections and DVs. We have shown that some of the moments used in the 
literature, such as the family of moments t/;^^'*^) [see eq. (4.2)], employed e.g. by ALEPH 
and OPAL [3, 5, 6], are perturbatively unstable and do not provide good approximations to 
the Borel resummed results of the models investigated. In the best of the cases, FOPT is 
capable of approaching the Borel results only at higher orders that are not available in the 
real QCD case. In addition, the ambiguity introduced by the power corrections is of the 
same order as the perturbative contribution, rendering the extraction of and condensates 
from these moments unreliable (similar observations were made in ref. [7]). Moments that 
contain an x term lie in between the two latter groups. In the reference model, where the 
contribution of the leading IR pole is significant, they display run-away behaviour in the 
perturbative series, which engenders larger uncertainties and potential instabilities in as 
results. The contribution from D = 4 in the OPE can also be quite sizeable depending on 
the coefficient of the x term in the weight function. This makes the use of these moments 
for determinations of as and condensates somewhat problematic, although the evidence 
against their use is less compelling than in the case of moments with only higher powers 
of x. 

Overall, regarding the FOPT/CIPT comparison, the moment analysis confirms the 
conclusions drawn from the inclusive tau hadronic width [18] in the following sense: when- 
ever perturbatively well-behaved moments are considered, FOPT shows better behaviour 
in the reference model that is believed to incorporate the main known features of large- 
order behaviour in QCD. CIPT underestimates the resummed value and therefore leads to 
systematically larger as values. 

In section 5, we have provided further evidence to the plausibility of the reference model 
suggested in ref. [18]. The matching procedure with the first two terms of a polynomial in 
u in eq. (3.6) was justified based on the behaviour of the Adler function on the complex 
plane and on comparisons with the large- /3o limit. This limit provides a laboratory for 
renormalon models, since the exact result to all orders is known. We have shown that 
models containing solely the leading singularities capture the general features of the exact 
result surprisingly accurately. What is more, the use of these models is sufficient to decide 
upon the best prescription for the RG improvement of the perturbative series. It therefore 
appears that the reference model is solid and survives criticisms raised in the literature. 

Acknowledgments 

MB is supported in part by the Gottfried Wilhelm Leibniz programme of the Deutsche 
Forschungsgemeinschaft (DFG). The work of DB is supported by the Alexander von Hum- 
boldt Foundation. MJ is partially supported by CICYTFEDER-FPA2008-01430, FPA2011- 



-32- 



25948, SGR2009-894, and the Spanish Consohder-Ingenio 2010 Program CPAN (CSD2007- 
00042). 



References 

[1] Y. Amhis et al. [Heavy Flavour Averaging Group], Averages of b-hadron, c-hadron, and 
tau-lepton properties as of early 2012, [arXiv:1207.1158 [liep-cx]]. 

[2] E. Braaten, S. Narison, and A. Pich, QCD analysis of the r hadronic widths 
Nucl. Phys. B 373 (1992) 581. 

[3] R. Barate et al. [ALEPH Collaboration], Measurement of the spectral functions of axial - 
vector hadronic tau decays and determination of as{m^), Eur. Phys. J. C 4 (1998) 409. 

[4] S. Schael et al. [ALEPH Collaboration], Branching ratios and spectral functions of tau 

decays: Final ALEPH measurements and physics implications, 
Phys. Kept. 421 (2005) 191 ]hep-ex/0506072]. 

[5] M. Davier, S. Descotes-Genon, A. Hocker, B. Malaescu, and Z. Zhang, The Determination of 
as from Tau Decays Revisited, 

Eur. Phys. J. C 56 (2008) 305 [arXiv:0803.0979 [hep-ph]]. 

[6] K. Ackerstaff et al. ]OPAL Collaboration], Measurement of the strong coupling constant 
and the vector and axial-vector spectral functions in hadronic tau decays, 
Eur. Phys. J. C 7 (1999) 571 ]hep-ex/9808019]. 

[7] K. Maltman and T. Yavin, ag{M^) from hadronic tau decays, 
Phys. Rev. D 78 (2008) 094020 [arXiv:0807.0650 [hep-ph]]. 

[8] D. Boito, O. Cata, M. Golterman, M. Jamin, K. Maltman, J. Osborne, and S. Peris, 

A new determination of ag from hadronic t decays, 

Phys. Rev. D 84 (2011) 113006 )arXiv:1110.1127 ]hep-ph]]. 

[9] D. Boito, M. Golterman, M. Jamin, A. Mahdavi, K. Maltman, J. Osborne, and S. Peris, 
An Updated determination of as from r decays, 
Phys. Rev. D 85 (2012) 093015 [arXiv: 1203.3146 [hep-ph]]. 

[10] B. Blok, M. A. Shifman and D. X. Zhang, 

An illustrative example of how quark-hadron duality might work, 

Phys. Rev. D 57 (1998) 2691 [Erratum-ibid. D 59 (1999) 019901] [hep-ph/9709333]. 

[11] I. I. Y. Bigi, M. A. Shifman, N. Uraltsev, A. I. Vainshtein, Heavy flavor decays, OPE and 

duality in two-dimensional 't Hooft model, Phys. Rev. D 59 (1999) 054011 [hep-ph/9805241]. 

[12] M. A. Shifman, Quark hadron duality, ]hep-ph/0009131]. 

[13] O. Cata, M. Golterman, S. Peris, Duality violations and spectral sum rules, 
JEEP 0508 (2005) 076 [hep-ph/0506004]. 

[14] O. Cata, M. Golterman, S. Peris, Unraveling duality violations in hadronic tau decays, 
Phys. Rev. D 77 (2008) 093006 [arXiv:0803.0246 [hep-ph]]. 

[15] O. Cata, M. Golterman, S. Peris, Possible duality violations in tau decay and their impact on 
the determination of as, Phys. Rev. D 79 (2009) 053002 larXiv:0812.2285 [hep-ph]]. 

[16] M. Jamin, What two models may teach us about duality violations in QCD, 
JHEP 1109 (2011) 141 [arXiv: 1103.2718 [hep-ph]]. 



-33- 



[17] M. Jamin, Contour-improved versus fixed- order perturbation theory in hadronic tau decays, 
JHEP 0509 (2005) 058 [hcp-ph/0509001]. 

[18] M. Beneke and M. Jamin, and the r hadronic width: fixed-order, contour-improved and 
higher-order perturbation theory, JHEP 0809 (2008) 044 [arXiv:0806.3156 [hep-ph]]. 

[19] A. A. Pivovarov, Renormalization group analysis of the tau-lepton decay within QCD, 
Z. Phys. C 53 (1992) 461 [Sov. J. Nucl. Phys. 54 (1991) 676] [Yad. Fiz. 54 (1991) 1114] 

[hep-ph/0302003]. 

[20] F. Le Diberder and A. Pich, Testing QCD with tau decays, Phys. Lett. B 289 (1992) 165. 

[21] P. A. Baikov, K. G. Chctyrkin and J. H. Kiihn, Order QCD Corrections to Z and r 
Decays, Phys. Rev. Lett. 101 (2008) 012002 ]arXiv:0801.1821 ]hep-ph]]. 

[22] S. Menke, On the determination of as from hadronic tau decays with contour-improved, fixed 
order and renormalon- chain perturbation theory, [arXiv:0904.1796 [hep-ph]]. 

[23] I. Caprini and J. Fischer, ag from r decays: contour-improved versus fixed-order summation 

in a new QCD perturbation expansion, 

Eur. Phys. J. C 64 (2009) 35 [arXiv:0906.5211 [hep-ph]]. 

[24] S. Dcscotes-Gcnon, B. Malaescu, A note on renormalon models for the determination of 
as(Mr), JarXiv: 1002.2968 Jhcp-ph]]. 

[25] G. Cvetic, M. Loewe, C. Martinez and C. Valenzuela, Modified Contour- Improved 
Perturbation Theory, Phys. Rev. D 82 (2010) 093007 [arXiv: 1005.4444 [hep-ph]]. 

[26] G. Abbas, B. Ananthanarayan and I. Caprini, Determination of as{M^) from Improved 

Fixed Order Perturbation Theory, Phys. Rev. D 85 (2012) 094018 [arXiv: 1202.2672 [hep-ph]]. 

[27] W. Marciano and A. Sirlin, Electroweak radiative corrections to t decay, 
Phys. Rev. Lett. 61 (1988) 1815. 

[28] E. Braatcn and C. S. Li, Electroweak radiative corrections to the semihadronic decay rate of 
the T lepton, Phys. Rev. D 42 (1990) 3888. 

[29] J. Erler, Electroweak radiative corrections to semileptonic tau decays. 
Rev. Mex. Fis. 50 (2004) 200 [hep-ph/0211345]. 

[30] I. S. Towner and J. C. Hardy, The evaluation of Vud and its impact on the unitarity of the 
Cabibbo-Kobayashi-Maskawa quark-mixing matrix, Rept. Prog. Phys. 73 (2010) 046301. 

[31] S. G. Gorishnii, A. L. Kataev, and S. A. Larin, The 0{a^) corrections to 

o'tot(e^e^ hadrons) and r(r~ i^r + hadrons) in QCD, Phys. Lett. B 259 (1991) 144. 

[32] L. R. Surguladze and M. A. Samuel, Total hadronic cross-section in e^e~ annihilation at the 
four-loop level of perturbative QCD, Phys. Rev. Lett. 66 (1991) 560. 

[33] M. Beneke, Renormalons, Phys. Rept. 317 (1999) 1 [hep-ph/ 9807443]. 

[34] P. Ball, M. Beneke, and V. M. Braun, Resummation of (/SqcCs)" corrections in QCD: 
Techniques and applications to the r hadronic width and the heavy quark pole mass, 
Nucl. Phys. B 452 (1995) 563 [hep-ph/9502300]. 

[35] M. Ncubcrt, QCD analysis of hadronic r decays revisited, 
Nucl. Phys. B 463 (1996) 511 [hep-ph/9509432]. 

[36] M. Beneke, Large order perturbation theory for a physical quantity, 
Nucl. Phys. B 405 (1993) 424. 



-34- 



[37] D. J. Broadhurst, Large-N expansion of QED: Asymptotic photon propagator and 

contributions to the muon anomaly, for any number of loops, Z. Phys. C 58 (1993) 339. 

[38] A. H. Mueller, On the Structure of Infrared Renormalons in Physical Processes at 
High-Energies, Nucl. Phys. B 250 (1985) 327. 

[39] M. Beneke, Renormalization scheme invariant large order perturbation theory and infrared 
renormalons in QCD, Phys. Lett. B 307 (1993) 154. 

[40] D. R. Boito, O. Cata, M. Golterman, M. Jamin, K. Maltman, J. Osborne, and S. Peris, 
Duality violations in tau hadronic spectral moments, 
Nucl. Phys. Proc. Suppl. 218 (2011) 104 [arXiv: 101 1.4426 [hep-ph]]. 

[41] A. Pich and J. Prades, Strange quark mass determination from Cabibbo-suppressed r decays, 
JHEP 10 (1999) 004 [hep-ph/9909244]. 

[42] E. Gamiz, M. Jamin, A. Pich, J. Prades, and F. Schwab, Determination of rUs and \ Vus\ 
from hadronic r decays, JHEP 01 (2003) 060 [hep-ph/0212230]. 

[43] E. Gamiz, M. Jamin, A. Pich. J. Prades. and F. Schwab, Vus and rUg from hadronic t decays, 
Phys. Rev. Lett. 94 (2005) 011803 ]hep-ph/0408044]. 

[44] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, QCD and resonance physics, 
Nucl. Phys. B 147 (1979) 385. 

[45] M. Jamin, Flavour- symmetry breaking of the quark condensate and chiral corrections to the 
Gell-Mann-Oakes-Renner relation, Phys. Lett. B 538 (2002) 71 [hep-ph/0201174]. 

[46] J. Beringer et al. (Particle Data Group), Phys. Rev. D 86 (2012) 010001. 



-35- 



